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Oscillating spin-orbit interaction in two-dimensional superlattices: 
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We consider ballistic transport through a lateral, two-dimensional superlattice with experimentally 
realizable, sinusoidally oscillating Rashba-type spin-orbit interaction. The periodic structure of the 
rectangular lattice produces a spin-dependent miniband structure for static SOI. Using Floquet 
theory, transmission peaks are shown to appear in the mini-bandgaps as a consequence of the 
additional, time-dependent SOI. A detailed analysis shows that this effect is due to the generation 
of harmonics of the driving frequency, via which e.g., resonances that cannot be excited in the case of 
static SOI become available. Additionally, the transmitted current shows space and time-dependent 
partial spin-polarization, in other words, polarization waves propagate through the superlattice. 

PACS numbers: 73.23.-b, 72.25.-b, 71.70.Ej 


I. INTRODUCTION 

Ballistic transport phenomena governed by time- 
dependent potentials are of fundamental interest mainly 
due to their close relation to important time-dependent 
quantum mechanical scattering effects. On the other 
hand, the possibility of controlling the electron dynam¬ 
ics using time-dependent gate voltages may result in 
practical applications. Additionally, as it has been 
demonstrated experimentally, spin-dependent properties 
of semiconductor materials - that are of exceptional im¬ 
portance e.g. in the context of spin-based quantum me¬ 
chanical information processingii^ - can also be con¬ 
trolled by gate voltagesi^ These results motivated us 
to investigate how oscillating Rashba-type spin-orbit in¬ 
teraction (SOI)^ affects spin-dependent conductance in 
two-dimensional super lattices. 

Spin-dependent transport phenomena in lateral super- 
lattices^^— have been investigated experimentally, mainly 
in a two-dimensional network of quantum ringsi^i^ Con¬ 
trol of spin geometric phase in semiconductor quantum 
rings has also been demonstratedi ^ Here, we focus on 
the ballistic regime and consider rectangular geometries, 
i.e., networks that consist of linear quantum wire seg¬ 
ments as building blocks (see Fig.[T|). The (quasi)periodic 
geometry of these devices results in a Rashba spin-orbit 
interaction controlled miniband structure)^ with char¬ 
acteristic energies orders of magnitude below the usual 
band widths. This is a direct consequence of the dif¬ 
ference between the lattice constant a (see Fig. [IJ, the 
order of 100 nm, and typical atomic separations. Since 
the position, width and even the existence of the non¬ 
conducting energy ranges (i.e., the mini-bandgaps) can 
be controlled experimentally via the strength of the SOI 
interaction, the conductance of the device is found to be 
tunable even at room temperature.— 

In the current paper we demonstrate that the time de¬ 
pendence of the spin-orbit interaction gives rise to new 
physical phenomena, leading to observable transmission 


peaks in the mini-bandgaps. We consider the combina¬ 
tion of oscillating and static SOI, where the latter one 
determines the miniband structure, while the oscillating 
part induce time-dependent effects. Note that transport 
related problems with oscillating SOI (but without the 
miniband structure) have been studied in Refs, (islll^ for 
a ring, and in Ref. [I^ for a ring-dot system. Application 
of Floquet’s theoryi^ allows us to obtain nonperturba- 
tive results, high order harmonics of the SOI oscillation 
frequency appear in the transmission. Floquet scatter¬ 
ing matrix theory is proven to be a useful mathematical 
tool for the description of periodically time dependent 
phenomena in diverse mesoscopic samples— Specifi¬ 
cally, several studies have discussed resonant phenomena 
of quantum dots and nanowires in the presence of time- 
dependent potential, see e.g. Refs [23l - [^ . 

Here, we will show that from a detailed analysis, we 
find the higher harmonics of the SOI oscillation frequency 
are responsible for the transmission peaks in the mini- 
bandgaps e.g., by allowing the excitation of resonances 
that are not coupled to the input/output leads for static 
SOI. 

The present paper is organized as follows. In section 
im we describe the model and methods that are used in 
the following. Physical consequences of this model are 
analyzed in Sec. HI. We present and discuss our time 
averaged results in section IV. Time resolved spin and 
charge density oscillations are discussed in Sec. V, while 
Sec. VI contains the summary and conclusions. 


II. MODEL AND METHODS 

The building blocks of the superlattices shown in Fig.[T] 
are linear, narrow (one dimensional) quantum wires. The 
corresponding Hamiltonian with Rashba-type SOI can be 


2 



for any two time instants r and r', allows us to calculate 
the time evolution for an arbitrary initial state. Con¬ 
cretely, the evolution operator, for which 

C/(fc, r) [3>(fc, s, r = 0)] = $(fc, s, r) (6) 

for any linear combination $(fc, s,t = 0) = Q;$i(fc,s) -I- 
/3d>2(A:,s), can be calculated explicitly: 


FIG. 1: Schematic view of the rectangular two-dimensional 
superlattice. Quantum wires represented by gray lines are 
subjected to oscillating Rashba-type SOI. Black arrows show 
the input and output leads in which no spin-orbit coupling is 
present. 


written^l^ as 


H{t) = h{l 
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n(a X e 
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4Q2 ’ 


. ( 1 ) 

where the unit vector n = cos + sin ipey points to 
the chosen positive direction along the wire, s denotes 
the dimensionless position variable (measured in units of 
the lattice constant, a) and we introduced the charac¬ 
teristic kinetic energy hil = It?/ 2m*a? . Note that anal¬ 
ogous Hamiltonian has also been used previously for 
quantum ringsi^^— The strength of the SOI is given by 
uj{t) = air)/a, where t = fit is the (dimensionless) time 
variable, and a denotes the Rashba parameter. The time 
dependence of a;(r) is assumed to be given by 


uj{t) = Wo -I- Wi cos(i^Q,r), (2) 


where Va is the frequency of the gate voltage oscillation 
(in units of Q). 


A. Floquet states and quasienergies 

Considering the solution of the time-dependent 
Schrddinger equation 


U{k,T)=e '^x 1 cos ( -I-sin(z/aT) ) 4- 


n 


ftVn 


. I ‘^0 J , Wife . 
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( 7 ) 


Here, 1 denotes the 2x2 identity matrix, while = 
cos LpGx + sin LfGy is the Pauli matrix corresponding to 
the direction of the actual lead, with ip representing the 
appropriate polar angle [for n = {By), <p = Q (7r/2)]. 

The time evolution operator © is diagonal in the basis 
of 

V>±(fc, s, r) = ^ ^ , (8) 


where 


e±(fc,T) = (/c^ ± ^fc)r =F sin(i^ar). (9) 

Note that the eigenspinors of the spin operator Oy, have 
the form 


l4>±) - ^2 ( ±*6*“^ ) ■ 

This means that the time-dependent basis spinors can be 
written as 

ip±{k,s,T) = U{k,T) [ip±ik,s,0)] = 

( 11 ) 

The exponential factor above can be factorized: 


,cl<i)(s, r) 
* dr 




( 3 ) 


it is seen that for an infinite, narrow quantum wire, any 
initial state can be expanded as a linear combination (a 
continuous one in case of the spatial variable) of spinor 
valued wave functions 


$i(fc,s) = 6*'=" 



$2(fc,s) = e*'=" 



( 4 ) 


which are expressed in the eigenbasis of For a given 
value of k (measured in units of 1/a), the action of the 
Hamiltonian on the states ([4]) becomes relatively simple, 
since the spatial derivatives have to be replaced by a 
multiplication by ik. This, together with the fact that 

( 5 ) 


g-ie±(fe,r) _ 


where 


( 12 ) 


el{k) = k^±'^k, (13) 

and the second term (that is periodic in time) can be 
expanded as: 

^ ^ (14) 

n— — oo ^ 

The Jacobi-Anger identity above (where Bessel functions 
of the first kind^ appear), explicitly shows that the 
states tjj±{k,s,T) can be called the Floquet states of the 
problem corresponding to the quasienergies of e±(fc). 


[i4(r),i4(r')]=0. 
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B. Global solution of the transport problem 


solving Eqs. (1201) for the wave number fc, we obtain that 


Having obtained the ’’time-dependent eigenspinors” 
((8|) of the Hamiltonian o, we have the solution of the 
time-dependent Schrddinger equation (|31) for an arbitrary 
initial state in an infinite quantum wire. The transport 
problem of the superlattice, however, involves quantum 
wire segments, and the solution has to obey appropriate 
boundary conditions. 

We assume that a monoenergetic plane wave input en¬ 
ters the network from the left hand side (see Fig. [1}: 

r) = (^ly (15) 

where, in dimensionless units, 

kin{E) = y/E. (16) 

The oscillating part of the SOI can induce high harmon¬ 
ics of frequency Va, leading to ’’sidebands” or Floquet 
channels^ in the transmission at dimensionless energies 


k'sliiEm) - ^ 1 / 4 ^"^^’ 

4o/(^-) = ^ ± (21) 

where the first two solutions correspond to the upper, 
while k^sQi to the lower sign in Fq. (l20ll . and the sub¬ 
script reminds us that these relations are valid in do¬ 
mains with oscillating SOI interaction. [Note the differ¬ 
ence between Fqs. (HI and (ED).] Combining Fqs. 

([20)) and ([Sp , we see that a general solution of the time- 
dependent Schrddinger equation ([3]) that involves the fre¬ 
quency components ED can be written in the following 
form: 

00 4 

'4’soi{s,t) = E E “Vi {ksOliEm), S,t) , (22) 

m— — oc i—1 

where the coefficients will have to be determined using 
the boundary conditions, and 


En=Eo+nVa, (17) 

with n being integer. Note that although SOI oscillations 
are obviously not quantized, this expression resembles the 
scenario when a quantum system absorbs/emits energy 
quanta proportional to i>a from/into a quantized field. 
The appearance of the frequencies ED in the time evo¬ 
lution of the quantum state of the system means that e.g., 
the reflected spinor valued wave function can be written 
as 


i '^+ {ksoA^rn),S,T) for 1=1,2 
s,r) = < 

[ ’>P-iksoii^rn),s,T) for f = 3,4. 

(23) 

Using the Jacobi-Anger expansion (11^1) . we obtain e.g.: 
i’m,! [k, S, r) =V'm,l {k, s, 0) 

1 — — 00 ^ ' 


'tprefl{s,T) 


^il-kir^iEr^)s-EnT)] 

n 


Tin \ 
r2n J ’ 


(18) 


where the coefficients Tin and r 2 n will be determined by 
the boundary conditions. Similarly, for the transmission 
(in the output arm) 


tptransis.T) = 


XE„)s-Er,r)] 


iln \ 

^2n J 


(19) 


As we shall see, the complete transport problem can be 
solved by imposing appropriate boundary conditions at 
the junctions - in frequency domain, that is, for each 
frequency component ED separately. By investigating 
Fqs. (I51)- (ITT)) . one can see that the time evolution of the 
solutions ([5]) involves a given frequency En, whenever 


el{k) = Em, (20) 


where m and n can be either equal or different. [In fact, 
once a term exp(—zUnx) appears in the time evolution 
of a state given by Eq. all other frequencies Em = 
En + (to — play a role - although it is possible that 
their weight in the Fourier expansion is negligible.] By 


At this stage, we obtained solutions to the time- 
dependent Schrddinger equation in all spatial domains: 
1pin{x,T) + lljrefl{x,T) iu the input arm, i:trans{x,T) in 
the output arm, and ipSOi{s,T) (with appropriate orien¬ 
tation, i.e., value of (p) inside the network. These spinor 
valued wave functions contain coefficients that are still 
to be determined. Using the coordinate system shown in 
Fig.m we require Re(/ci„) > 0, Im(fci„) < 0 in the input 
arm, and Re(fci„) > 0, Im(fci„) > 0 for the transmit¬ 
ted solution, in order to ensure left propagating reflected 
waves, right propagating transmitted waves and evanes¬ 
cent solutions that decay as a function of the distance 
from the central region. Since the domains on which the 
functions ipsoiiSjT) are defined are finite, and propaga¬ 
tion in both directions is possible, there are no restriction 
for kgQj- At the junctions (input, output and internal 
ones) we apply Griffith’s boundary conditions^ for each 
frequency component separately. In this way the spinor 
valued wave function is continuous at any time instant at 
all the junctions, and the net spin current density that 
leaves/enters any given junction disappear always. The 
resulting infinite system of linear equations (for more de¬ 
tails, see the Appendix) can be truncated. Practically, for 
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the results presented in this paper, it turns out that con¬ 
sidering approximately 50 values of (n = —25,... 25) 
is sufficient to achieve accurate results. This can be 
checked reliably via calculating the time averaged trans¬ 
mission and reflection probabilities (see the next section) 
that has to add up to unity. When this sum is not close 
enough to 1 (within 10“^ relative error), we increase the 
number of frequency components that are taken into ac¬ 
count. 



Eo/hn 

FIG. 2: (Color online) Conductance (in units of Go = ^|-) 
of a 7 X 7 array in the presence of oscillating and stationary 
Rashba-type spin-orbit interaction [see Eq. ([2j]. Panel a) fo¬ 
cuses on an energy range that corresponds to a mini-bandgap 
without the oscillating part of the SOI (see the inset). In 
this panel, G is plotted for different oscillating SOI ampli¬ 
tudes oJi, with ojo being fixed. According to panel b), when 
ojo/fl is considerably lower than oJi/H, the mini-bandgap dis¬ 
appears. Parameters are: uJo/Q. = 1.0, Va = 3.0 (panel a)) 
and uji/Q, = 0.4, Va = 3.0 (panel b)). 


III. OBSERVABLES 


the solution 


V’(s,t) 


Mi(s,t) \ 
M 2 (s,r) J 


(25) 


[That stands for tpin{x,T) + lfjrefl{x,T), tptrans{x,T) Or 

i’SOiis^r), depending on the location]. The position de¬ 
pendent (unnormalized) electron density is given by 


n(s,r) = tp\s,T)ijj{s,T) 


|wi(S:'r)P + |w2(s,r)p, (26) 


while the density for the spin-up and spin-down electrons 
(in the z direction) reads 

n-f(s,T) = |ui(s,t)P, n|(s,T) = |M2(s,r)p. (27) 

Focusing on the spin degree of freedom, one can construct 
the quantum mechanical spin density operator 

1 f |ui(s,t)P Mi(s,r)u^(s,r) \ 

^ ’ n(s,T) V'«i('S,t)m2(s,t) |u2(s,t)P y’ 

(28) 

which is defined only for nonzero electron density. Note 
that Tr[/9] = 1 (by construction), and for spin polarized 
states Tr[p^] = 1 also holds. However, when we would like 
to perform calculations for completely unpolarized input, 
the easiest way is to consider two different ipin states sep¬ 
arately, with their spinor parts being antiparallel, and fi¬ 
nally add the results incoherently, with equal statistical 
weight (i.e., 1/2) being associated to each states. For¬ 
mally, this is equivalent to an input spin density operator 
that is 1/2 times the 2x2 identity matrix [see Eq. (IMl) ]. 
In this case Tr[pf^] = 1/2, suggesting that the quantity 
Tr[p^(s,T)] can be an appropriate local measure of spin- 
polarization (As it is often used in different contexts, see 
e.g. Ref. [4l|). However, in our case it is more intuitive 
to express the degree of spin polarization as the length 
of the vector 


/ Tr[p{s,T)a^]\ 

S(s,r) = Tr[p(s,r)cry] (29) 

V T'r[p{s,T)a^] J 

that describes the spin orientation (as the components 
are the expectation values of the Pauli matrices). As it 
can be shown easily, S is a unit vector for spin polarized 
states, while its length is zero when p is proportional to 
unity. The ’’purity” 


p = VM (30) 

will be used to measure the degree of spin polarization. 
(Note that p G [0,1].) 

As an important quantity that does not depend on 
time, we calculate the time averaged transmission prob¬ 
ability. 


r 


1 

/ 'Joutix 

^in Jq 


0, t)dt, 


The results to be presented in the current section are 
related to physical quantities that can be calculated using 


■T 


(31) 
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where the usual quantum mechanical probability current 
densities appear, and T = 27r/z/Q,. Explicit calculation 
using Eqs. (USD and (TT^ shows that the time averaged 
conductance, which is proportional to T, can be written 
in units of 2e^ jh as: 

^ ~ T\ 12 I \hl2\l. ( p 'I + l^2n| )kin{En)- 

{\a\^ + \b\^)hn{Eo) ^ 

(32) 



FIG. 3: (Color online) Conductance (in units of Go = ^|-) 
of a 7 X 7 array as a function of the input energy Eq and the 
frequency Va of the oscillating part of the SOI. Additional 
parameters: cuo/fl = 1.0, cuo/H = 0.3. 


IV. TIME AVERAGED CONDUCTANCE 
PROPERTIES IN THE MINI-BANDGAPS 

Fig. [2] shows the time averaged conductance [see 
Eq. (IS2D] S'® 3- function of the input energy Eq for dif¬ 
ferent parameters. The role of the oscillating part of the 
SOI is clearly seen by comparing the inset and the main 
plot in panel (a): in the energy range corresponding to a 
mini-bandgap of the system with wi = 0, G is practically 
zero for the static system, but the oscillating part of the 
SOI induces several conductance peaks. In the following 
we will analyse these peaks. 

Let us recalli^ that there are no bandgaps without 
static SOI (wq = 0), and - as a rule of thumb - the 
widths of these energy ranges with zero conductance in¬ 
crease as a function of [wqI. As panel b) of Fig. [5] shows, 
sufficiently strong oscillating SOI can smear out the mini¬ 
band structure, even in cases when mini-bandgaps would 
still exist for wi = 0. This is related to the fact that the 
larger the magnitude of wi is, the more pronounced the 
peaks in Fig. [2] (a) are: broader and higher peaks in the 
mini-bandgaps lead to the disappearance of the bandgap 
itself. As we shall see later, the peaks seen in Fig. [2] (a) 
are related to the conductance via the Floquet channels 
corresponding to the harmonics (HZD, i.e., their appear¬ 
ance is a nonlinear effect. Therefore they are expected to 
play a more important role as the amplitude of the SOI 
oscillations - that generate them - increases. 


The position and physical origin of the conductance 
peaks in the mini-bandgaps needs a more detailed anal¬ 
ysis. To this end it is instructive to see the dependence 
of the position of the peaks on the frequency Va of the 
oscillating SOI. Fig. [3] shows G{Eo) for different values 
of Va- Intuitively, based on their widths, heights, and 
shapes, we can identify three kinds of peaks in the mini- 
bandgap shown in Fig. [3] Local maxima that are simi¬ 
lar in this sense are plotted using the same symbols and 
colors in this figure. As we shall see later, visual sim¬ 
ilarity corresponds to similar physical interpretation as 
well. First, focusing on the projections on the bottom 
plane, we can see that the position of the local conduc¬ 
tance maxima (M) changes linearly with z/q,. More con¬ 
cretely, M{va) = const-|-nz/ct, where n has integer values 
that differ only in sign for the peaks that are plotted 
using the same symbols in Fig. [3l Although the ’’driv¬ 
ing field” (Rashba-type SOI) is completely classical (the 
oscillations are not assumed to be quantized), this lin¬ 
earity, together with Fig. SI resembles the process of 
emission/absorption of oscillation quanta by the quan¬ 
tum system - i.e., the electron. Consequently, Figs. 0] 
(a), (c), (d) and (f) - that shows the weight of the fre¬ 
quency components CZD in the output give by Eq. dUD 
~ can be interpreted relatively easily. For panel (a) [(c)] 
the maximum of the transmission is shifted by (twice) 
Vet below Eq. That is, the emergence of harmonics (side¬ 
bands around the input energy Eq) allows ’’mapping” of 
energy ranges with nonzero conductance into the mini- 
bandgap. Moreover, the position of peak (a) [(c)] has 
an energy (note that we are using dimensionless units) 
distance of Va {2va) from the lower band edge around 
Eo/hn = 150 (see Fig. |3D. Similarly, peaks (d) and 
(f) ’’map” the nonzero conductance (that can be seen 
in Fig.|2](a) above Eq/MI= 165) by the ’’absorption” of 
one [(f)] or two [(d)] ’’oscillation quanta” Va- 

The narrow peaks whose energy distributions are de¬ 
noted by (b) and (e) in Fig. 0] are unrelated to the band 
edges, and their dependence on the driving frequency is 
different from that of peaks (a), (c), (d) and (f). Accord¬ 
ing to Fig. 01 the energy difference between these peaks is 
2va , to a very good approximation (within 2% relative er¬ 
ror). As we have checked by singular value decomposition 
(SVD) of the matrix that describes the boundary condi¬ 
tions (fitting at the junctions) for constant SOI (wi = 0), 
the energy value in the middle of these two peaks corre¬ 
sponds to a strong, multiply degenerate singular value. 
In other words, there are solutions that can be added 
to the scattering problem, that is, the global spinor val¬ 
ued wave function is not uniquely determined. However, 
these singular solutions have the property, that the cor¬ 
responding electron densities are zero at the input junc¬ 
tion. (E.g., for the parameters corresponding to Fig. 0] 
the singularity appears at Eq/KQ. = 157.66, and there are 
standing probability waves around the input junction -- 
with a node being at this point - so that the character¬ 
istic wavelength is a/2.) This means that these solutions 
are ’’closed”, have no coupling to the input lead. In other 
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FIG. 4: (Color online) The weights + |timp of the fre¬ 

quency components Em [see Eqs. dia and m] in the trans¬ 
mitted spinor valued wave function as a function of the har¬ 
monic order (the integer m). Figures (a),(b) •. • (f) corre¬ 
spond to the 6 peaks that can be identihed in the front curve 
{va = 3) of Fig. [3] - from left (lowest energy) to the right 
(highest energy), respectively. For the clear identification of 
the peaks, the same symbols were used as in Fig. [3] Addi¬ 
tional parameters are; wo/fl = 1.0, ujq/Q, = 0.3. 


words, singular solutions cannot be excited directly, i.e., 
not at the energy value where the matrix describing the 
boundary conditions is indeed singular without the os¬ 
cillating part of the SOI. That is why the weight of the 
corresponding frequency component is practically zero in 
Figs, m (b) and (e). However, when nonlinear effects in¬ 
duced by the oscillating part of the SOI give rise to high 
harmonics, the corresponding wavelengths do not all re¬ 
sult in destructive interference at the input junction: the 
conductance becomes nonzero. 

Thus, all the peaks that appear in the mini-bandgap 
are related to the emergence of high harmonics of the 
frequency of the driving SOI oscillations, but the de¬ 
tailed physical mechanisms are different for the broad 
and narrow local conductance maxima. In the first case 
the edges of the mini conduction bands are ’’mapped” 
into the mini-bandgap, while a strong, narrow resonance 
is being excited in the latter case. 


V. CHARGE DENSITY OSCILLATIONS AND 
SPIN POLARIZATION 


Time-resolved details of the transmission can be vi¬ 
sualized by plotting snapshots of the electron density 
along the network for various time instants. Here we 
consider completely unpolarized input, i.e., a plane wave 


with completely random spin polarization: 


Pin{x,T) 


}i„i[ki„{Eo)x-EoT)] 

2 

}_Akir,{Eo)x-Eor)] 1 0 ^ 

r [oij- 



0 0 Y 

0 

(33) 


As we can see in Fig. [5l the density n{x, y, r) has sev¬ 
eral maxima around the input junction, the location of 
which oscillates periodically during a cycle of duration 
T = 27r/r'a. This figure reveals an additional difference 
between the broad and the narrow peaks that appear in 
the mini-bandgaps: in the latter case [e.g., panels (b) and 
(e) in Fig. 0], the excitation of an internal resonance of 
the network results in considerably higher particle densi¬ 
ties. 





FIG. 5: (Golor online) Snapshots of the electron density along 
a 7 X 7 array for time instants indicated in the labels of the 
vertical axes. Note that n is not normalized. [As a reference: 
the value of n = 1 corresponds to the input plane wave (in this 
case with unpolarized spin state).] Panels (a)-(e) correspond 
to the first, broad peak seen in Fig. [3] [panel (a) of Fig. [4], 
while panel (e) correspond to the parameters of the second, 
narrow peak in Fig. [3] [panel (b) of Fig. [4]. 


Considering the output, it is instructive to point out 
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that SOI oscillations can lead to temporal spinpolariza- 
tion. This effect, in simpler geometries without mini¬ 
band structure, has already been demonstrated,— thus 
our current findings - besides providing a more detailed 
physical interpretation of the effect - indicate that this 
is a general consequence of the time-dependent SOI, be¬ 
ing essentially independent from the geometry of the sys¬ 
tem. Additionally, let us emphasize, that temporal spin- 
polarization is closely related to the oscillation of the SOI, 
no polarization appears for the case of static SOI. With¬ 
out the time dependence of the spin-related properties 
of the device, strong, symmetry-based considerations^ 
related to the equilibrium spin currents rule out spin po¬ 
larization effects. 

Figure [5] demonstrates that in the output arm, the de¬ 
gree of spin polarization characterized by p can be close 
to unity such that the electron density is still nonzero. 
Let us note that spin polarization and density fluctua¬ 
tions appearing in this figure propagate away from the 
network in a wave-like manner. The arrows in Fig. |6] (b) 
represent the spin orientation (l2^ separately for the two, 
opposite input spin direction [see the first line of Eq. (IMl) ] 
the incoherent sum of which constitutes the input density 
matrix (1331) . More precisely, the arrows visualize the spin 
direction in a local coordinate system, they point from 
(x, 0, 0) to (x + Sx, Sy, Sz)- By investigating both panels 
of this figure, one can see the physical origin of the po¬ 
larization effect: the spin directions corresponding to the 
two different input spinors rotate in a different way, they 
are not always antiparallel (which is the case for static 
SOI). In fact, there are spacetime points when these di¬ 
rections are almost the same, resulting in a remarkable 
partial polarization, p. This emphasises the role of the 
oscillating part of the SOI in the spin polarization effect 
shown in Fig. |6l 


VI. CONCLUSIONS 

We developed a model for the description of time and 
spin-dependent transport phenomena in rectangular, lat¬ 
eral superlattices. Motivated by recent experimental 
possibilities, the combined effect of static and oscillat¬ 
ing Rashba-type spin-orbit interaction (SOI) were con¬ 
sidered. The static part of the coupling induces an ex¬ 
perimentally controllable miniband structure, while the 
oscillating part gives rise to transmission peaks in the 
mini-bandgaps. This effect is general for networks that 
contain more than 5x5 junctions. In order to see a 
clear miniband structure with relatively low computa¬ 
tional costs - without loss of generality - we have chosen 
the size of 7 x 7 to demonstrate our results. We identified 
the physical mechanisms responsible for the appearance 
of conductance peaks in the mini-bandgaps, and have 
shown that the heights and positions of these peaks can 
be controlled by the amplitude and frequency of the SOI 
oscillations. These observations may lead to e.g. narrow 
band, controllable energy filters. 



FIG. 6: (Color online) Panel a): The degree of spin polariza¬ 
tion p and the particle density n in the output arm of a 7 x 7 
array at r = 0. The spin orientation corresponding to ’’spin- 
up” and ’’spin-down” (in the z direction) input spinors are 
shown in panel b). Parameters are: uJo/Q = 1.0, ojo/fl = 0.3, 
Eo/Ml = 153.5. 
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Appendix 

As an example, let us consider the output junction 
(where the output lead is connected to the network). 
Using So to denote the location of this point, Griffith’ 
boundary conditions^ require the solution to be contin¬ 
uous at sq. That is, all the neighboring spinor valued 
wave functions evaluated at this point should be equal 
at any time instant. As an example, considering a quan¬ 
tum wire segment that joins the output junction and the 
output wire itself, we can write: 

^SOl{^ — — '4^trans{^ — 

X] { So, t )) 

m 

+a'^'tp+{kgoi{Em), So, r))] 

+ So,t)) 

-haV-(fcso/(^m), So, t) )] I 

_ ^gi[fcin(Br.)so-BT.'r)] f tln \ ^ 

rj ' 






























where the probability amplitudes a*, the wavenumber 
kgoii^m) and the states '0±(A:gQj(£'m), Sq, r)) are in¬ 
troduced in Eqs. m, (Eu, m and (12311 ■ 

The periodicity (in time) of the problem offers a rela¬ 
tively simple way to imply the condition above, since (via 
Fourier series expansion) it is possible to work in the fre¬ 
quency domain. E.g., for the frequency component Em 
we have: 


l<^+) 

I 

yitZ/o; J y ^2m\-^m) J 

(34) 


The second part of the boundary conditions is related 


to the quantum mechanical probability current density, 
which, in the presence of the SOI, reads: 


J(s,r) = 2Re 


. d uj{t) 




(35) 


where \E'(s,t) denotes a solution to the time-dependent 
Schrodinger equation ([3]). (A derivation that leads to an 
analogous expression for a ring, can be found in Ref. [i^.l 
As one can check, having continuity imposed [Eqs. (IMl) ]. 
the condition that the net current density that flows in a 
junction (or, depending on the sign, out of it) should be 
zero at any time instant;^ turns into a set of linear equa¬ 
tions involving spatial derivatives. With an appropriate 
truncation of the infinite system of equations describing 
the boundary conditions, a global solution of the scatter¬ 
ing problem can be achieved. 
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